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ABSTRACT 

The recently suggested swing interaction be- 
tween fast magnetoso ni c an d Alfven waves 
l)Zaaarashvili fc RobertsL |2002a) is generalized 
to inhomogeneous media. We show that the fast 
magnetosonic waves propagating across an applied 
non-uniform magnetic field can parametrically 
amplify the Alfven waves propagating along the field 
through the periodical variation of the Alfven speed. 
The resonant Alfven waves have half the frequency 
and the perpendicular velocity polarization of the 
fast waves. The wavelengths of the resonant waves 
have different values across the magnetic field, due 
to the inhomogeneity in the Alfven speed. There- 
fore, if the medium is bounded along the magnetic 
field, then the harmonics of the Alfven waves, which 
satisfy the condition for onset of a standing pattern, 
have stronger growth rates. In these regions the 
fast magnetosonic waves can be strongly 'absorbed', 
their energy going in transversal Alfven waves. We 
refer to this phenomenon as 1 Swing Absorption" . 
This mechanism can be of importance in various 
astrophysical situations. 

Key words: MHD waves; wave coupling. 
1. INTRODUCTION 

Wave motions play an important role in many 
astrophysical phenomena. Magnetohydrodynamic 
(MHD) waves may transport momentum and energy, 
resulting in heating and acceleration of an ambient 
plasma. A variety of waves have recently been de- 
tected in the solar atmosphere using the SOHO and 
TRACE spacecraft. Hence, an understanding of the 
basic physical mechanisms of excitation, damping 
and the interaction between the different kinds of 
MHD wave modes is of increasing interest ijRobertsl 
l2f)04j) . Formally speaking, there is a group of di- 
rect mechanisms of wave excitation by external forces 
(e.g. turbulent convection, explosive events in stel- 
lar atmospheres, etc.) and wave dissipation due 
to non-adiabatic processes in a medium (such as 
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viscosity, thermal conduction, magnetic resistivity, 
etc.). There is a separate group of wave amplifica- 
tion an d damping processes due to re sonant mech- 
anisms j|GoossensL 11991: Pocdts, 2002J). This means 
that particular wave modes may be damped (ampli- 
fied) due to energy transfer (extraction) into or from 
other kinds of oscillatory motions, even when wave 
dissipation is excluded from consideration. 

Recently, a new kind of interaction between dif- 
ferent MHD wave modes, bas ed on a pa r amet- 
ric action, has been suggested llZaaarashviil 12001 : 
I Zaaarashvili fe Robertsl)2002alpl : IZaaarashviri et al. . 
2002, 2004). In this case the mechanism of wave 
interaction originates from a basic physical phe- 
nomenon known in classical mechanics as 'paramet- 
ric resonance', occurring when an external force (or 
oscillation) amplifies the oscillation through a pe- 
riodical variation of the system's parameters. The 
mechanical analogy of this phenomenon is a mathe- 
matical pendulum with periodically varying length. 
When the frequency of the length variation is twice 
the frequency of the pendulum oscillation, the ampli- 
tude of the oscillation grows exponentially in time. 
Such a mechanical system can consist of a pendulum 
(transversal oscillations) with a spring (longitudinal 
oscillations). A detailed de scription of such a me- 
chanica l system is given in [Zaqarashvili fc Roberts! 
(2002^) (hereinafter referred to as Paper I). 

Here we present a brief report on our recent results 
on swing interactions of fast and Alfven modes in 
inhomogeneous media (Shergelashvili et al. (2004)). 



2. BASIC EQUATIONS AND EQUILIB- 
RIUM MODEL 

Consider a magnetized medium with zero viscosity 
and infinite conductivity, where processes are as- 
sumed to be adiabatic. Then the macroscopic dy- 
namical behaviour of this medium is governed by the 
ideal magnetohydrodynamic (MHD) equations: 

?E + v.( p u) = o, (l) 
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Figure 1. Schematic view of the inhomogeneous back- 
ground configuration and the directions of wave prop- 
agation (polarization). 
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where p and p are the plasma pressure and density, 
U is the velocity, B is the magnetic field of strength 
B = \B\ and 7 denotes the ratio of specific heats. 

We consider an equilibrium magnetic field directed 
along the z axis of a Cartesian coordinate system, 
Bq = (0, 0, Bq(x)). The equilibrium magnetic field 

Bq and density po are inhomogeneous in x. The force 
balance condition, from Eq. J2Jl, gives the total (ther- 
mal + magnetic) pressure in the equilibrium to be a 
constant: 



Po(x) 
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= constant. 



Also, the equation of state 

p {x) =p {T (x),po(x)), 



(5) 



(6) 



relates the unperturbed pressure, density po(x) and 
temperature Tq(x). 

Eqs. (l)-(4) are linearized around the static (C/o = 0) 
equilibrium state (5). This enables the study of the 
linear dynamics of magnetosonic and Alfven waves. 
A schematic view of the equilibrium configuration is 
shown in Figure 1. 

3. FAST MAGNETOSONIC WAVES 

Now let us consider the propagation of "pure' fast 
magnetosonic waves across the magnetic flux sur- 
faces, i.e. along the x axis, taking u = (u x ,0,0) and 



d/dy = 0, d/dz = 0. The equation governing the 
dynamics of fast magnetosonic waves propagating 
across the magnetic field lines in the inhomogeneous 
medium can be obta ined from equations JU-® as 
flRobertsllT98Hll99ll) : 
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47r J dx 
where u x is the velocity perturbation 



u 2 po(x)u(x) = 0. 



(7) 



The solution of this equation for different particular 
equilibrium conditions can be obtained either ana- 
lytically or numerically. Equation Q describes the 
propagation of a fast wave with speed (C 2 + V^) 1 ^ 2 , 
where C s — ("/Po/po) 1 ^ 2 is the sound speed and 
V A = {Bfi/Airpo) 1 / 2 is the Alfven speed. 

We study Eq. Q in accordance with the boundary 
conditions: 

u(0) = u{L x ) = 0, (8) 

corresponding to fast waves bounded by walls located 
at x — and x — L x . These boundary conditions 
make the spectrum of fast modes discrete: Eq. (JJJ 
has a nontrivial solution only for a discrete set of 
frequencies, 

w = w n =w ,wi, (9) 

In this case the solutions for different physical quan- 
tities can be represented as 

u x = av(x) sin(w„i), p\ = ar(x) cos(w n i), 

b z = ah(x) cos(w n i), (10) 

where the density and magnetic field perturbations 
are related to the velocity perturbations through 



1 / dp dv(x) 
r( , ) = _^ (x) _ +po __ 

h(x) = ± Lx) d -^ + B ?p 
u!„ v dx dx 
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Here (and subsequently) we use a subscript n to de- 
note the frequency of a given standing fast mode. 
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Now consider Alfven waves that are linearly polar- 
ized in the y direction and propagate along the mag- 
netic field (see Figure 1). In the linear limit these 
waves are decoupled from the magnetosonic waves 
and the equation governing their dynamics is: 

d 2 b 

" (13) 



dt 2 
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where Va(x) is Alfven speed. It is clear from Eq. (|13J) 
that the phase speed of this mode depends on x para- 
metrically. Therefore, an Alfven wave with a given 
wave length propagates with a "local' characteristic 
frequency. Each magnetic flux surface can evolve in- 
dependently in this perturbation mode. 



4.1 



Propagating Alfven waves 



Let us now address the non-linear action of the fast 
magnetosonic waves, considered in the previous sec- 
tion, on Alfven waves. We study the weakly non- 
linear regime. This means that the amplitudes of the 
fast magnetosonic waves are considered to be large 
enough to produce significant variations of the en- 
vironment parameters, which can be felt by prop- 
agating Alfven modes, but too small to affect the 
Alfven modes themselves. Hence, the magnetic flux 
surfaces can still evolve independently. Therefore, 
as in paper I, the non-linear terms in the equa- 
tions arising from the advective derivatives u x db y / dx 
and (po + Pi) u x duy/dx are assumed to be negligi- 
ble. Under these circumstances the governing set of 
equations takes the form (see Paper I): 
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These equations describe the parametric influence of 
fast magnetosonic waves propagating across the mag- 
netic field on Alfven waves propagating along the 
field. An analytical solution of Eqs. ifl^jl and (fT5|) is 
possible for a standing pattern of fast magnetosonic 
waves, the medium being assumed bounded in the x 
direction. 

Combining Eqs. (|14J) and Ijl5|l we obtain the follow- 
ing second order partial differential equation: 
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Writing 



by = h y (z, t) exp 

we obtain 
d 2 h y 1 
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Finally, applying a Fourier analysis with respect to 
the z coordinate, 



hy(z,t) 



hy{k Z: t)e z dk Z: 



(22) 



and neglecting the second and higher order terms in 
a, we obtain the following Mathieu-type equation: 
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dt 2 
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It should be noted that the expression l|20|l for 
S2(x,t) consists only of terms of second order and 
higher in a, and so can be neglected directly for the 
case of weakly non-linear action addressed here. 

Equation H23|l has a resonant solution when the fre- 
quency of the Alfven mode to a is half of uj n : 



lu a = k z V A (x) w -uj t , 



(25) 



This solution can be expressed as 



h y (k z ,t) = h y (k z ,t = 0)es 
where 



cos — t — sin — t 
2 2 



5(x) = akiVi(x)F(x). 
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The solution has a resonant nature within the fre- 
quency interval 



(28) 



Similar expressions have been obtained in the Pa- 
per I for a homogeneous medium. In that case, the 
Alfven speed is constant and, therefore, the fast mag- 
netosonic waves amplify the Alfven waves with the 
same wavelength everywhere. In the case of an in- 
homogeneous Alfven speed, the resonance condition 
(I25|l implies that the wavelength of the resonant har- 
monics of the Alfven waves depends on x. This 
means that the fast magnetosonic waves now am- 
plify Alfven waves with different wavelengths (but 
with the same frequency) in different magnetic flux 
surfaces (i.e., different a;- values). 

4.2 Standing Alfven waves 

When we consider a system that is bounded in the z 
direction, the boundary conditions along the z axis 
introduce an additional quantization of the wave pa- 
rameters. In particular, in this case each spatial har- 
monic of the Alfven mode can be represented as 



hy{km> t) COs(k m z) , 



(29) 



where k m = nm/L z (m = 1,2,...) and L z is the 
characteristic length of the system in the z direction. 
This then leads to a further localization of the spatial 
region where the swing transfer of wave energy from 
longitudinal to transversal oscillations is permitted. 
The resonant condition (|25|l implies that 



n in V A (x n, m J 
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Therefore, the resonant areas are concentrated 
around the points x n ^ m for which condition (|30|l 
holds. Within these resonant areas the longitudi- 
nal oscillations damp effectively and their energy 
is transferred to transversal oscillations with wave 
numbers k z = k m satisfying the resonant conditions. 
These resonant areas are localized in space and can 
be referred to as regions of swing absorption of the 
fast magnetosonic oscillations of the system. The 
particular feature of this process is that the energy 
transfer of fast magnetosonic waves to Alfven waves 
occurs at half the frequency of the fast waves. 
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Table 1. Values of the constant parameters used in 
the calculation of our illustrative solutions. The di- 
mension of C is g 1 ! 2 cm^^l 2 ' 1 /dyn 1 ^ 1 



For the configuration described by the equilibrium 
profiles 1)3 l|l - l)33|l the resonant condition (|30|l yields 
the areas of swing absorption (located along the x 
axis) as the solutions of the following equation: 



u n C poo +Poi 
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In Figure (panels Ai and A2) we plot 



(35) 



5. NUMERICAL SIMULATION 

In this section we consider in detail the process of 
swing absorption of fast waves into Alfven waves. 
We consider a numerical study of equation sub- 
ject to the boundary conditions JSJ. We study as 
an example, the case of a polytropic plasma when 
both the thermal and magnetic pressures are linear 
functions of the x coordinate: 
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where, poo, Poi> ^-ooi ^01 an d C are constants, and 
L x denotes the length of the system along the x di- 
rection. The pressure balance condition (j5j) immedi- 
ately yields 

Poi = —z-- (34) 

The solution of the wave equation depends on the 
values of the above set of constant parameters. In 
general, different equilibrium regimes can be consid- 
ered including those corresponding to different limits 
of the plasma [3: /3 -C 1, [3 w 1 and (3 > 1. The 
values of all constant parameters are given in Ta- 
ble n We took arbitrary values of parameters, but 
they are somewhat appropriate to the magnetically 
dominated solar atmosphere (say the chromospheric 
network). In Figure we show the profiles of av(x) 
for the standing wave solutions, for two cases with 
different modal 'wavelength'. Panels A and B, re- 
spectively, correspond to the characteristic frequen- 
cies: ujo ~ 1-87 ■ 10~ 2 s" 1 (period 5.61 min) and 
UJ2 ~ 5.66 • 10~ 2 s^ 1 (period 1.85 min). 
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(solid line) and 
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(dotted line) against the normalized coordinate 
x/L x . These curves correspond to the zeroth-order 
harmonic of the fast magnetosonic mode shown in 
Figure (panel A) and the standing Alfven mode 
with wave numbers m ~ 3 (panel Ai) and m = 4 
(panel A 2 ). 

In order to examine the validity of the approxima- 
tions we made during the analysis of the govern- 
ing equation (|18)l . we performed a direct numeri- 
cal solution of the set of equations (|14fl - i|15|l and 
obtained the following results. The Alfven mode 
m = 3 is amplified effectively close to the resonant 
point x/L x « xq.3/L x — 0.7967. This is shown 
on panel Ci of Figure |3| Far from this resonant 
point, the swing interaction is weaker, as at the point 
x/L x w xoa/ L x — 0.0767 (panel Bi). For the Alfven 
mode m = 4 we have the opposite picture: the area of 
'swing absorption' is situated around the point £0,4 
(see panel B2, Figure 0) and the rate of interaction 
between modes decreases far from this area, as at 
point xo,3 (panel C2, Figure|!]J). In these calculations 
we took a — 0.015. 

Similar results are obtained for the fast magnetosonic 
mode shown in panel B of FigureEl corresponding to 
a = 0.03. In this case, the fast magnetosonic mode 
effectively amplifies four different spatial harmonics 
of Alfven modes, viz. m = 8, 10, 11 and 12 (for details 
see Shergelashvili et al. 2004). A similar analysis 
can be performed for the case of any other equilib- 
rium configuration and corresponding harmonics of 
the standing fast magnetosonic modes. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

x/L x/L 

X X 



Figure 2. Sample solutions of the standing fast magnetosonic modes. Panel A: zeroth-order harmonic n = 0, 
period T = 5.6069 min., a = 0.015 ; Panel B: second-order harmonic n = 2, period T = 1.8514 min., a = 0.03 




Figure 3. The swing absorption of the fast mode n = by the standing Alfven modes m = 3 and 4- 



6. DISCUSSION AND CONCLUSION 

The most important characteristic of swing absorp- 
tion is that the velocity polarization of the amplified 
Alfven wave is strictly perpendicular to the veloc- 
ity polarization (and propagation direction) of fast 
magnetosonic waves. This is due to the parametric 
nature of the interaction. For comparison, the well- 
known resonant absorption of a fast magnetosonic 
wave can take place only when it does not propagate 
strictly perpendicular to the magnetic flux surfaces 
and the plane of the Alfven wave polarization. In 
other words, the energy in fast magnetosonic waves 
propagating strictly perpendicular (i.e. k± = 0) to 
the magnetic flux surfaces cannot be resonantly 'ab- 
sorbed' by Alfven waves with the same frequency 
polarized in the perpendicular plane. This is be- 
cause the mechanism of resonant absorption is anal- 
ogous to the mechanical pendulum undergoing the 
direct action of an external periodic force. This force 
may resonantly amplify only those oscillations that 
at least partly lie in the plane of force. On the con- 
trary, the external periodic force acting parametri- 
cally on the pendulum length may amplify the pen- 
dulum oscillation in any plane. A similar process 
occurs when the fast magnetosonic wave propagates 
across the unperturbed magnetic field. It causes a 
periodical variation of the local Alfven speed and 
thus affects the propagation properties of the Alfven 
waves. As a result, those particular harmonics of 
the Alfven waves that satisfy the resonant conditions 
grow exponentially in time. These resonant harmon- 
ics are polarized perpendicular to the fast magne- 
tosonic waves and have half the frequency of these 
waves. Hence, for standing fast magnetosonic waves 
with frequency uj n , the resonant Alfven waves have 
frequency ~u; n /2. 

In a homogeneous medium all resonant harmonics 
have the same wavelengths (see Paper I). Therefore, 
once a given harmonic of the fast and Alfven modes 
satisfies the appropriate resonant conditions (Eqs. 
(23) and (25) in Paper I), then these conditions are 
met within the entire medium. Thus, in a homo- 
geneous medium the region where fast modes effec- 
tively interact with the corresponding Alfven waves 
is not localized, but instead covers the entire sys- 
tem. However, when the equilibrium is inhomoge- 
neous across the applied magnetic field, the wave- 
lengths of the resonant harmonics depend on the 
local Alfven speed. When the medium is bounded 
along the unperturbed magnetic field (i.e. along the z 
axis), the resonant harmonics of the standing Alfven 
waves (whose wavelengths satisfy condition (|30|l for 
the onset of a standing pattern) will have stronger 
growth rates. This means that the 'absorption' of 
fast waves will be stronger at particular locations 
across the magnetic field. In the previous section we 
showed numerical solutions of standing fast magne- 
tosonic modes for a polytropic equilibrium (j>q ~ p^) 
in which the thermal pressure and magnetic pres- 
sure are linear functions of x. Further, we performed 



a numerical simulation of the energy transfer from 
fast magnetosonic waves into Alfven waves at the 
resonant locations, i.e. the regions of swing absorp- 
tion. 

The mechanism of swing absorption can be of impor- 
tance in a variety of astrophysical situations. 
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